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£N| ■ In this work, we have discussed the Maxwell's electrodynamics in non-linear forms in FRW uni- 

verse. The energy density and pressure for non-linear electrodynamics have been written in the 
electro-magnetic universe. The Einstein's field equations for flat FRW model in loop quantum cos- 
mology have been considered if the universe is filled with the matter and electro-magnetic field. We 
separately assumed the magnetic universe and electric universe. The interaction between matter 
r—{ _ and magnetic field have been considered in one section and for some particular form of interaction 

term, we have found the solutions of magnetic field and the energy density of matter. We have 
also considered the interaction between the matter and electric field and another form of interaction 
1 term has been chosen to solve the field equations. The validity of generalized second law of thermo- 

^Nj . dynamics has been investigated on apparent and event horizons using Gibb's law and the first law 

of thermodynamics for magnetic and electric universe separately. 
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I. INTRODUCTION 

• i— i ■ 

Observations of the redshift of supernovae type la [l], 0] and cosmic microwave background 0, 0] show 
that our universe is expanding with acceleration and lead to the search for a new type of matter which 
violates the strong energy condition, known as dark energy. The dark energy has the property that the 
energy component produce sufficient negative pressure, which drives the cosmic acceleration. There are many 
candidates supporting this behavior. Cosmological constant A is the most popular candidate of dark energy 
satisfying the equation of state (EoS) pa = —PA- There are other strong favored candidates like quintessence, 
which is composed by a scalar field [f| Q with self-interacting potential. The EoS of the fluid distribution 
of the universe is given by p = up, where w is known as the EoS parameter. The universe will accelerate if 
u) < —1/3. u> = —1 represents ACDM model and uj < — 1 corresponds to the phantom dominated model (tJ. 
Recently, many candidates with variable EoS play the crucial role of dark energy to drive the acceleration of 
the universe namely, Chaplygin gas d, [t| , Tachyonic field [l(| , holographic dark energy [ll| , agegraphic dark 
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energy [ljj], Ricci dark energy [l3j . Hessence [13], DBI-essence [lj|, K-essence [ll|, dilaton dark energy []jj etc. 

At present, we live in an epoch where the dark energy and the dark matter are comparable. For this 
purpose, the interacting dark energy models have been studied to explain the cosmic coincidence problem 
18]. Till now, many works have been proposed for this interacting dark energy In recent years, the 

model of interacting dark energy has been explored in the framework of loop quantum cosmology (LQC). 
The LQC is the application in the cosmological context of loop quantum gravity (LQG) [2l| - [23| . which is 
a theory trying to quantize the gravity with a non-perturbative and background independent method. By 
studying the early universe inflation and the fate of future singularity in LQC, it is found that the big bang 
singularity, the big rip singularity and other future singularities can be avoided [2J]. It has been verified that, 
the cosmological evolution in LQC for quintessence model is same as that in classical Einstein cosmology, 
whereas for the phantom dark energy, the loop quantum effect significantly reduce the parameter spacetime 
required by stability. Recently, the dynamics of phantom, quintom and hessence dark energy models in LQC 
have been studied [25j ]. 

On the other hand, a new approach [2t| has recently been taken to avoid the cosmic singularity through 
a non-linear extension of the Maxwell's electromagnetic theory. The associated Lagrangian and the resulting 
electrodynamics can theoretically be justified, based on different arguments. The homogeneous and isotropic 
non-singular FRW solutions can be obtained [26] by considering a generalized model of Maxwell's electrody- 
namics, where local covariant and gauge-invariant Lagrangian depend on the field invariants up to the second 
order, as a source of classical Einstein's equations. Exact solutions of the Einstein's field equations coupled to 
non-linear electrodynamics (NLED) may hint at the relevance of the non-linear effects in strong gravitational 
and magnetic fields. An inhomogeneous and anisotropic non-singular model for the early universe filled with 



Born- Infold type non-linear electromagnetic field was studied 27]. Recently, there are several works on NLED 
in various situations 12814311 . 



In Einstein's gravity, the connection between black hole thermodynamics and Einstein's equations was first 
discovered in [32J by deriving the Einstein's equations from the proportionality of entropy and horizon area 
together with the first law of thermodynamics. The thermodynamical laws also have been applied in the 
cosmological context, considering universe as a thermodynamical system bounded by the apparent horizon. At 
the apparent horizon, the first law of thermodynamics is equivalent to the Friedmann equations and generalized 
second law (GSL) is obeyed at the horizon. There are several studies in thermodynamics for dark energy filled 
universe on apparent and event horizons [331- 37]. Ther mody namical properties in non-linear electrodynamics 
has extensively studied in [38j]. In our previous work [39j . we have assumed the apparent, event, Hubble 
and particle horizons of the FRW universe and have studied the first law and GSL of thermodynamics in 
non-linear electrodynamics with magnetic field only. In the present work, we extend the previous work in LQC 
by considering both the electric and magnetic fields and examine the validity of GSL in apparent and event 
horizons of the magnetic universe and electric universe. 

The investigation has been done in the following way : We have briefly discussed the Maxwell's electrodynam- 
ics in non-linear forms in section II. The energy density and pressure for non-linear electrodynamics have been 
written in the electro-magnetic universe. The Einstein's field equations for flat FRW model in loop quantum 
cosmology have been considered if the universe is filled with the matter and electro-magnetic field. We sepa- 
rately assumed the magnetic universe (E = 0) and electric universe (B = 0). The interaction between matter 
and magnetic field have been considered in section III and for some particular form of interaction term, we have 
found the solutions of magnetic field and the energy density of matter. In section IV, we have considered the 
interaction between the matter and electric field only and another form of interaction term has been chosen 
to solve the field equations. In section V, the validity of generalized second law of thermodynamics have been 
investigated on apparent and event horizons using Gibbs' law and the first law of thermodynamics for magnetic 
universe and electric universe separately. Finally, some conclusions are drawn in section IV. 

II. BASIC EQUATIONS IN NON-LINEAR ELECTRODYNAMICS 

The Lagrangian density in Maxwell's electrodynamics (linear) can be written as [3l| 



(1) 
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where F^ v is the electromagnetic field strength tensor and /xo is the magnetic permeability. The generalization 
of Maxwell's electro-magnetic Lagrangian (non-linear) up to the second order terms of the fields is given by [3l[ 

C = — F + loF 2 + r/F* 2 (2) 

where lj and rj are arbitrary constants and 



F* = F* V F^ U (3) 

where F* v is the dual of F^. Here we consider the FRW model of our universe. Since the spatial section of 
FRW geometry are isotropic, electromagnetic field can generate such a universe only if an averaging procedure is 
performed (40ll4l|. In this situation, the energy density and the pressure of the NLED field should be evaluated 
by averaging over volume. So the corresponding energy-momentum tensor for non-linear electro-magnetic theory 
has the form 



The modified Lagrangian in non-linear electrodynamics for accelerated universe is considered as [42j 



C = -- F + aF 2 +/3F- 1 (5) 

where a and /3 are arbitrary (constant) parameters. As seen this Lagrangian contains both positive and 
negative powers of F. The second (quadratic) term dominates during very early epochs of the cosmic dynamics, 
while the Maxwell term (first term above) dominates in the radiation era. The last term is responsible for 
the accelerated phase of the cosmic evolution [43[ ■ The above Lagrangian density yields a unified scenario to 
describe both the acceleration of the universe (for weak fields) and the avoidance of the initial singularity, as a 
consequence of its properties in the strong-field regime. 

The energy density and pressure for electro-magnetic (EM) field are given by [42| 



PF 



-C~AE 2 C F 



(G) 



and 



p F = C--{2B 2 ~E 2 )C F 



(7) 



where, B and E are respectively magnetic field and electric field. Now, the electro-magnetic field has the 
expression F = 2(i3 2 — E 2 ), so the explicit forms of the energy density and the pressure for electro- magnetic 
field will be (H H 



PF = ^B 2 +E 2 ) 



4a{B 2 - E 2 )(B 2 +3E 2 ) 



P(B 2 - 3E 2 ) 
2(B 2 - E 2 ) 2 



(8) 



and 



VF = + E 2 ) -^(B 2 - E 2 )(5B 2 - E 2 ) + f*:%$ 0) 

Due to the loop quantum effect, the standard Friedmann equation in LQC can be modified by adding a 
correction term as (for fiat model) [4^, |46| 
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8?T G Ptotai 
ti — — - — Ptotai 1 

Pi 



(10) 



and 



H = -47T G(ptotal + Ptotai) ( 1 



^Ptotai 
Pi 



(11) 



Also the energy-conservation equation is given by 



Ptotai + 3H(p tota l + Ptotai) = 



(12) 



Pm + Pf where, p m and p m are energy density and pressure for matter 
iPm, Pi = V3n 2 -f 3 G 2 h is called the critical loop quantum density, 7 is 



with, p^tai = p m + Pf and p to tal 
obeying the equation of state p m = w, 
the dimensionless Barbero-Immirzi parameter and H = | is the Hubble parameter. 

Now if we consider the homogeneous electric field E in plasma, it gives rise to an electric current of charged 
particles and then rapidly decays. So the squared magnetic field B 2 dominates over E 2 , i.e., in this case, the 
average value < E 2 and hence F — 2B 2 . So F is now only the function of magnetic field (vanishing 
electric component) and hence the FRW universe may be called the magnetic universe. Similarly, if the 
average value < B 2 >m 0, then F = —2E 2 . So F is now only the function of electric field (vanishing magnetic 
component) and hence the FRW universe may be called the electric universe. In the following two sections, we 
shall assume the magnetic and electric universes separately. 



III. INTERACTION BETWEEN MATTER AND MAGNETIC FIELD 

In the magnetic universe, we assume E = 0. Therefore the energy density and pressure are (using equations 
(8) and (9)) 



9B = \B 2 -^-^ 



(13) 



PB -Q B + ZB 2 



and subsequently 



(14) 



Ptotai = Pm + \B 2 - iaB 4 - 



2B 2 



(15) 



and 



1 1 r2 20a 4 7/3 

Ptotai — U m p m + -13 —tS + ^2 



(16) 



Now if we consider the interaction between matter and magnetic field, then the conservation equation (12) 
becomes 



and 



p m + 3H(1 + w m )p m = Q 



(17) 



5 



f>B + 3H(p B + PB ) = -Q (18) 
where Q is the interaction term. Now let us take the interaction component Q as 



Q = 2SB (l - WaB 2 + -j^j H 



(19) 



This expression is used to simplify the calculation procedure. Here 5 may be treated as interaction parameter 
which is a small positive quantity. Using the expressions of p B , Pb from (15) and (16) and the relation between 
p m and p m , the above two equations (17) and (18) give the solutions of B and p m in the following forms: 

B = -S + ^ (20) 

and 



p m = p a-^+^ + ™ 



a 6 {(3 + 5 4 - 16aS e ) ^ 16aB% 48a 2 aSB 2 



35 3 (1 + U m ) 3(1 - U m ) 3cj m - 1 



+ , f-3/3 + <5 4 - 48a«5 6 + 6(3 2 i^(l + 3c m ), 1, + A 

a z b A (l + 3u} m ) \ 2 2 B Q 

-4/3 a *k[i(l + 3^), 2, + fg] + ^ 2Fl[ l (1 + 3^), 3, ( 21 ) 
where -Bo and po are positive integration constants and 2-F1 is the hypergeometric function. 



IV. INTERACTION BETWEEN MATTER AND ELECTRIC FIELD 

This section deals with the case B = i.e., the universe is filled with matter and electric field. Thus the 
expressions of the energy density and pressure for the electric field become (using equations (8) and (9)) 

PE = \E 2 + l2aE i + ^L (22) 



and therefore 



PE ~6 ~ ~3~ ~ ~§EA (23) 



Ptotal = Pra + ^E 2 + \2aE 4 + ^ ( 24 ) 



and 



Ptotal = u mPm + ^E 2 - ^-E* - ^jL (25) 

In this case, we consider the interaction between matter and electric field. So the conservation equation (12) 
takes the form 



() 



p m + 3H(1 + w m )p m = Q (26) 

and 

PE + ZH{p E +p E ) = -Q (27) 
To simplify of the calculation, the interaction component Q is taken as 

Q = 35E 2 (l + l6aE 2 + -^j H (28) 

where 5 is a small positive quantity. In this case, exact solutions for p m and E from the above two equations 
can only be obtained if they are solved numerically using the expressions of pe and pe together with the 
relation between p m and p m . 



V. GENERALIZED SECOND LAW OF THERMODYNAMICS 

In this section, the validity of the generalized second law of thermodynamics is studied. It states that, the 
sum of entropy of total matter enclosed by the horizon and the entropy of the horizon does not decrease with 
time. In the following, we consider apparent and event horizons only. The variation of entropy inside the 
horizon will be calculated via Gibb's equation and the variation of entropy on the horizon will be calculated 
using first law of thermodynamics. Hence we shall examine the validity of GSL of thermodynamics of the 
magnetic universe and electric universe separately bounded by the above mentioned horizons. 



A. Apparent Horizon 

We know that radius of apparent horizon in the FRW universe, 

Ra = 4 (29) 



Therefore, 



H 2 

Case-I: E = 0: Using equations (10), (11), (15) and (16) 



H 



Ra = -4 (30) 



B 2 (l - IQaB 2 + -Jr) + §(1 + w m )p 
p m + ^~ 4aB 4 - 2g 

Considering the net amount of energy crossing through the apparent horizon in time dt as [4 



Ra ^ — ■b 1 _'_, (3D 



-dE = ^R\B{ P total + Ptotal)dt (32) 

and assuming the validity of first law of thermodynamics on the apparent horizon, i.e., 



— dE — T^dSA 



(33) 
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S A +S, 




Fig.l represents rate of change of total entropy of apparent horizon i.e., Sa + Si against the scale factor a for 
B = with interaction for w m ~ 1/3 (solid line), w m = (dotted line) and w m — —0.5 (dashed line). 



where, Sa and Ta are the entropy and temperature on the apparent horizon, we have 



dS A 
dt 

Again from the Gibb's eqn 



4ttR\H 
T a 



2B 2 



1 - 16aB 2 + 4i 
B 4 



+ (1 + w m )p„ 



(34) 



T A dSi = dEj + ptotaidV (35) 
where, Sj is the internal entropy, V is the volume and Ex = pV is the internal energy, we have 



dSj 
dt 



AirR 2 A 
Ta 



IB 1 



1 - \<auB l 



(R A - HR A ) 



(36) 



From eqns (34) and (36), the rate of change of the total entropy becomes 



dt 



(S A + Si) = 



AttR 2 a 
T a 



2B 2 



1 - IQaB 2 



B 4 



+ (1 + W m )pr, 



Ra 



(37) 



Substituting the expression (31) in eqn (37) and using the expressions of B and p m from eqns (20) and (21) 
respectively, we plot the rate of change of total entropy of the apparent horizon, i.e, Sa + Si against the scale 
factor in figure 1 with interaction (5 = 0.0001) for different matter components i.e., w m = 1/3 (solid line), 
w m = (dotted line) and w m = —0.5 (dashed line). From the figure, we see that the rate of change of total 
entropy for apparent horizon is initially positive but as a increases, it decreases and in late time, it becomes 
negative. So we may conclude that, under the loop quantum cosmological effects, the GSL is valid initially, but 
after certain stage of the evolution of the universe, the GSL will not valid for apparent horizon for interacting 
scenarios in the magnetic universe. 



Case-II : B = 0: 



Proceeding in the same way as in the previous case, the expression for the rate of change of the total entropy 
becomes 
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Fig.2 represents rate of change of total entropy of apparent horizon i.e., Sa + Si against the scale factor a for 
E = with interaction for w m — 1/3 (solid line), w m = (dotted line) and w m — —0.5 (dashed line). 



dt 



(S A + Si) 



4irR 2 A 



2E 2 



1 + WaE 2 



+ (1 + w m )p n 



Ra 



(38) 



Substituting the expression of Ra from eqn (30) and using the numerical solutions of E and p m , we plot the 
rate of change of total entropy of the apparent horizon, i.e, Sa + Si against the scale factor in figure 2 with 
interaction (S — 0.0001) for different matter components i.e., w m = 1/3 (solid line), w m = (dotted line) and 
w m = —0.5 (dashed line). From the figure, we see that the rate of change of total entropy for apparent horizon 
is always positive. Thus it may be concluded that even under the loop quantum cosmological effects, the GSL 
is always satisfied in late time for apparent horizon for interacting scenarios in the electric universe. 



B. Event Horizon 



The event horizon radius is given by 



R E = a 



da 



(39) 



The differential eqn of which can be written as 



Re — HRe — 1 

Considering the net amount of energy crossing through the event horizon in time dt as 



(40) 



- dE= 4xR%H(ptotal +Ptotai)dt 
and assuming the validity of first law of thermodynamics on the event horizon, i.e, 



(41) 



— dE = TsdSs 



(42) 
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Case-I : E=0: Fig. 3 represents rate of change of total entropy of event horizon i.e., Se + Si against the 
scale factor a for B — with interaction for w m = 1/3 (solid line), w m = (dotted line) and w m — —0.5 (dashed line). 



we have the rate of change of the total entropy as in the following cases: 



Case-I : E = 0: 



d ,„ „ . AttR% 



2B 2 f 2 8 

1 - 16<xB 2 + —j 

3 V B 4 



+ (1 + W m )pr, 



{HR E - 1) 



(43) 



Case-II : B = 0: 



1 + 16aE z 



^4 ) + 



(FEe - 1) 



(44) 



Substituting the expressions of H, B, E and p m in eqns (43) and (44), the rate of change of total 
entropy of the event horizon, i.e, Se + Si is plotted against the scale factor in figures 3 and 4 for the above two 
cases with interaction (S = 0.0001) for different matter components i.e., w m = 1/3 (solid line), w m — (dotted 
line) and w m = —0.5 (dashed line). From the figures, we see that the rate of change of total entropy for event 
horizon is always positive for both magnetic and electric universes and hence the GSL is always satisfied for 
event horizon for interacting scenarios of the magnetic and electric universes. 



VI. DISCUSSIONS 



In this work, we have briefly discussed Maxwell's electrodynamics in non-linear forms for accelerating 
universe. The energy density and pressure for non-linear electrodynamics have been written in magnetic 
universe in one section and for universe with electric field only in another section. The Einstein's field equations 
for loop quantum cosmological model have been considered for FRW model of the universe. The interaction 
between matter and electric and magnetic fields have been incorporated separately and for particular forms 
of interaction terms, we have found the solutions for both electric as well as magnetic fields and the energy 
density of matter. 
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Case-II : B=0 Fig. 4 represents rate of change of total entropy of event horizon i.e., Se + Si against the 
scale factor a for E = with interaction for w m = 1/3 (solid line), w m — (dotted line) and w m = —0.5 (dashed line). 



In addition to this, our endeavor was to investigate the validity of the generalized second law of ther- 
modynamics of the universe bounded by the apparent and event horizons. The variation of entropy 
has been calculated inside the horizon using Gibb's equation and that on the horizon using the first law of 
thermodynamics. After that, we have studied the GSL of the universe bounded by the above mentioned horizons. 

In figures 1-4, the variation of total entropy on the apparent and event horizons have been drawn 
against the scale factor a for interacting scenarios (5 = 0.0001) of magnetic universe as well as universe 
with electric field only, for w m = 0,1/3,-0.5. From figure 1, we see that the rate of change of total 
entropy was initially positive but in late epoch, it becomes negative and thus the GSL is not valid for 
the magnetic universe on the apparent horizon. Whereas on the event horizon for the magnetic universe, 
the GSL remains always valid. This can be clearly seen from figure 3. Figure 2 represents the rate of 
change of total entropy on the apparent horizon for the universe with electric field only, which again shows 
the validity of the GSL throughout the evolution of the universe for all the cases. Finally, from figure 4, 
we see that GSL was initially not valid on the event horizon, but in late time it is satisfied for all types of matter. 
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